Spin Hall effect at interfaces between HgTe/CdTe quantum wells and metals 
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We study the spin-dependent transmission through interfaces between a HgTe/CdTe quantum 
well (QW) and a metal for both the normal metal and the superconducting cases. Interestingly, 
we discover a type of spin Hall effect at these interfaces that happens to exist even in the absence 
of structure and bulk inversion asymmetry within each subsystem (i.e., the QW and the metal). 
Thus, this is a pure boundary spin Hall effect which can be directly related to the existence of 
exponentially localized edge states at the interface. We demonstrate how this effect can be measured 
and functionalized for an all-electric spin injection into normal metal leads. 

PACS numbers; 73.23.-b,73.63.-b,74.45.-|-c 
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I. INTRODUCTION 

The Spin Hall effect (SHE) is a rich physical phe- 
nomenon which, since its prediction has led to many 
interesting discoveries in the field of spintronics. Real- 
ized in nonmagnetic systems, the SHE allows for an all- 
electrical manipulation of spin. The underlying interac- 
tion of the SHE is spin-orbit coupling with the observable 
consequence that a transverse spin current is generated if 
an electrical charge current is driven in the longitudinal 
direction. Notably, this can happen due to impurity scat- 
tering [llli the extrinsic SHE, or due to band structure 
effects [3, Ij], the called intrinsic SHE. The latter case 
is usually directly related to either structure inversion 
or bulk inversion asymmetry of the underlying nanos- 
tructure, resulting in Rashba or Dresselhaus spin-orbit 
coupling, respectively. 

Remarkably, by now, both cases have been experi- 
mentally observed: the extrinsic SHE in semiconduc- 
tor heterostructures by optical Q and electronic [6] 
measurement techniques and the intrinsic SHE in two- 
dimensional semiconductors {?!] as well as in HgTe/CdTe 
heterostructures by combining the SHE with the so-called 
quantum spin Hall effect (QSHE) in a single device 
The QSHE is yet another type of spin Hall effect that 
exists at the boundary of a two-dimensional topological 
insulator realized in HgTe/CdTe quantum wells (QWs) 
[l0| . It refers to the existence of protected metallic edge 
states propagating in opposite directions and forming a 
single set of Kramers partners at each edge flll - [l3l | . The 
topological phase transition behind this effect is related 
to an inversion of bands with opposite parities due to 
the strong spin-orbit interaction and it is well described 
by a massive Dirac equation [l^ . The critical point is 
then reached when the mass of the Dirac fermions goes 
to zero. Directly at that point, the system behaves like a 
single valley Dirac fermion; this has been experimentally 
confirmed [14| . 

Recently, Yokoyama and co-workers predicted a giant 



spin rotation at a junction between a normal metal and 
QSHE system [l^; this is one of the main motivations of 
our work. Here, we go substantially beyond this predic- 
tion in two different ways. First and most importantly, 
we discover a different type of interface SHE at junc- 
tions between HgTe/CdTe QWs and metals (for both 
the normal metal and the superconducting cases). Sec- 
ond, we functionalize our findings to propose a device for 
all-electric spin injection into normal metal leads in the 
absence of ferromagnetic contacts. All our predictions 
apply to nanostructures in the ballistic transport regime. 

We show below that this type of SHE is intimately re- 
lated to the coexistence of propagating and evanescent 
modes at the interface between a QSHE system and a 
metal. Mathematically, this comes from the fact that 
the underlying low-energy Hamiltonian contains terms 
that are linear and quadratic in the electron wave vec- 
tor. Physically, this interplay can happen near a band- 
inversion crossing that in the case of HgTe/CdTe QWs 
drives the topological phase transition jl2l | . Interestingly, 
this effect exists even in the total absence of structure 
and bulk inversion asymmetry within each subsystem. 
To the best of our knowledge, this is the first prediction 
of a SHE in a composite system that does not break these 
symmetries in either subsystem. 
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FIG. 1: Junction between a HgTe/CdTe QW (light gray) and 
a metal (dark gray). Due to the applied voltage V , electrons 
are injected from a contact with an angle of incidence 6 with 
respect to the interface. At the interface they are either re- 
flected or transmitted. 

The paper is organized as follows. In Sec. II, we 
present the model Hamiltonian for a HgTe/CdTe quan- 
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turn well and describe the formalism needed to compute 
the scattering probabilities of a HgTe QW/metal junc- 
tion. In Sec. Ill, we investigate the spin- and angular- 
dependence of transmission as well as the Andreev reflec- 
tion probabilities and analyze their asymmetric behavior 
as an interfacial spin Hall effect. We propose to mea- 
sure this effect and functionalize it by an experimental 
realization described in Sec. IV. We conclude in Sec. V. 
Details of calculations are given in the Appendix. 



II. MODEL 

This section is devoted to the theoretical framework of 
the junction between a HgTe/CdTe quantum well and 
a metal. We present the model Hamiltonian for the 
HgTe/CdTe quantum well and describe the scattering 
matrix method which allows us to calculate to the scat- 
tering probabilities through the interface. 



A. Hamiltonian 

The band structure of a HgTe/CdTe QW is derived 
from the eight-band Kane model [16(1 . Near the topo- 
logical phase transition, it can be described by an ef- 
fective four-band Hamiltonian with two subbands com- 
monly called El and HI that have opposite parities. The 
El and the HI subbands are both doubly degenerate due 
to time- reversal symmetry (TRS). We refer to these de- 
generate Kramers partners as spin t, i- Near the T point, 
the effective Hamiltonian can be written as UM 



_ fh{k) 

^'^'^)- y h*{-k) J ' 



(1) 



where /i(k) = e(k)/2x2 + rfa(k)cr° represents the spin-t 
block in the {\E1 t) , \H1 1)) space, h*{-k) the spin-J, 
block in the (|£'l |) , \H1 |)) space, and cr° denote the 
Pauh matrices. In Eq. ([T|), e(k) = C - DP, da{k) = 
[Afc^, -Aky,M{k)], and M(k) = M - BP, where k = 
{kx,ky) is the in-plane momentum, P = + ky and 
M represents the Dirac mass. The parameters A, B, C, 
D, and M depend on the geometry of the HgTe/CdTe 
QW. The sign of M relative to B distinguishes the trivial 
phase (M > 0) from the non-trivial phase (M < 0, B 
being negative). The eigenenergies of ^(k) are given by 
E± = e(k) ± d{k) where d(k) = y/A^P + M(k)2 and ± 
refers to the conduction (valence) band. 



B. Scattering method 



X < and C = Cr for a; > 0. We make the reasonable 
assumption to model the normal metal as a highly doped 
HgTe/CdTe QW (similar to the treatment of the corre- 
sponding problem in graphene 17[) and the proximity- 
induced superconductivity with a step-like varying pair 
potential [lg| . We assume that the only boundary of the 
problem is the QW/metal interface. Therefore, we ana- 
lyze bulk state transport in the remainder of this paper. 

The scattering states in the bulk of the HgTe/CdTe 
QW can be written for each block of the Hamiltonian 
([T]) separately. They are plane wave (two-component) 
spinor wave functions that depend on all the parameters 
of the Hamiltonian, in particular on k^ and ky. Since 
we assume translation invariance in the y direction, the 
transverse wave vector ky is conserved. Due to the block 
structure of the Hamiltonian, we can analyze the spin-f 
scattering problem separately from the spin- J, case. In 
the QW region {x < 0), the spin- up quasiparticles are 
described by two-components spinors [15, 19, 20] 



^ f ±d(k) + Af (k) 

\ A[kx iky) 



(2) 



where the ± sign labels the conduction (valence) sub- 
band. The dispersion relation E{k) yields two possible 
wave vectors, namely 



fcl,2 = 



2(^2 -D2) 



1/2 



(3) 



with 7 = -A2 + 2MB - 2D{Cl - E) and t] = 4(5^ - 
£)2)[m2 — (Cl — E^]. We make the further assumptions 
that ^2 > £)2 and M < which are moti- 

vated by typical parameters for HgTe/CdTe QWs. Under 
this choice of parameters, we find that fci is real and k2 
is imaginary. This results in the coexistence of propagat- 
ing and evanescent modes on both sides of the boundary, 
which gives rise to the interface SHE further discussed 
below. On the left-hand side, propagating electrons have 
a real wave vector k^ — ±/ci cos 9, where 9 is the angle 
of incidence and ± labels the incident (reflected) mode. 
Meanwhile, evanescent electrons are described by a com- 
plex wave vector k^ = —iiky — fc|)~^/2^ 

At the metal side of the interface (a; > 0), the outgoing 
scattering states have a similar form as the ones for x < 
but the definitions of the wave vectors fci^2 contain the 
parameter Cr instead of Cl (see Appendix 1 for more de- 
tails). To obtain the scattering amplitudes for transmis- 
sion through the junction, we match the wave functions 
and their derivatives at a; = 0. Once this calculation has 
been done for the spin-f block the corresponding spin-J, 
problem follows by TRS. 



We consider a junction between a HgTe/CdTe QW and 
a normal metal or an s-wave superconductor as depicted 
in Fig. [T] The interface is assumed to be perfect and 
located at a; = 0. The model of this junction relies on a 
step-like variation of the bands, modeled by C = Cl for 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section, we investigate the scattering probabili- 
ties, which depend on the mass M and angle of incidence 
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6. We first focus on the junction between a HgTe QW 
and a normal metal for two different doping configura- 
tions. Afterward, we consider the contact of a HgTe QW 
and an s-wave superconductor involving Andreev reflec- 
tion processes. In both cases, the angular asymmetric 
behavior of the scattering probabilities is connected to 
the presence of evanescent modes. The coexistence of 
evanescent and propagating waves is necessary for the 
appearance of an interfacial spin Hall effect. 



A. Normal metal case 

We first consider the scattering amplitudes at the 
Fermi energy E — through a HgTe/CdTe QW-normal 
metal junction when both subsystems are n-doped. The 
transmission probability and the amplitude squared of 
the evanescent mode on the right-hand side of the junc- 
tion are shown in Fig. [5] as functions of the mass param- 
eter M and the angle of incidence 0. 




FIG. 2: (Color online) Angle dependence of the transmis- 
sion probability (left panel) and the amplitude squared of the 
evanescent mode on the metal side (right panel) for the spin-f 
block for a junction between a low-n-doped HgTe/CdTe QW 
and an n'-doped metal at the Fermi level [E — 0) as func- 
tions of the mass parameter M, where —Cl ~ 1 meV and 
— Cr = 500 meV. The other parameters are those found in 
typical HgTe QW experiments, i.e., A = 4 eV A, B = -70 
eV A^ and D = —50 eV A^. A clear asymmetry is present 
for negative values of M. 

For positive values of the gap M, the transmission 
probability exhibits a symmetric behavior with respect 
to 9 and has a rather weak amplitude. On the contrary, 
for a large negative mass M, this amplitude increases and 
tends to peak at negative 9. Comparing the transmission 
with the amplitude squared of the evanescent mode (see 
the right panel of Fig. [5]), a correlation between the two 
is visible. Hence, the magnitude of the amplitude of the 
evanescent mode is connected to the angle and spin de- 
pendence of the transmission. (Note that the spin- J, case 
just follows from the substitution 9 —9 in the results 
shown in Fig. O ) This gives rise to a spin current at 
the interface in the transverse direction [21|. Only the 
coexistence of propagating and evanescent modes at the 
interface gives rise to this effect. In the parameter regime 
where evanescent modes are absent, one can easily show 
that the transmission of the two spin blocks has to be 



symmetric with respect to 9. We remark in passing that 
the p — n' case yields similar results to the n — n' case 
discussed here. 




FIG. 3: (Color online) Angle dependence of the transmis- 
sion probability (left panel) and the amplitude squared of the 
evanescent mode on the metal side (right panel) for the spin-f 
block for a junction between a low-n-doped HgTe/CdTe QW 
and a p-doped metal at the Fermi level {E = 0) as functions of 
the mass parameter M, where —Cl~1 meV and Cr — 500 
meV. The values of the parameters A, B, and D are the same 
as in Fig. (2] Now, a clear asymmetry is present for values of 
M close to the critical point. 

Interestingly, when the metal side is p-doped (and not 
n'-doped as discussed before in Fig. [2]), a similar asym- 
metry in the spin- and angle- dependent transmission 
holds, but the maximum then appears near the critical 
point M ^ (see Fig. 13]). 



B. Superconducting case 

We now turn to the analysis of a junction between 
a HgTe QW and a superconductor (SC). Then, the 
Hamiltonian of Eq. ([IJ must be completed by particle- 
hole symmetry and the pairing potential matrix A(k). 
[We assume s-wave singlet pairing; see Eq. ^H)) in the 
Appendix.] In our model, the HgTe/CdTe QW on the 
left-hand side of the interface (x < 0) with a chemical 
potential Cl has no pairing potential, while the Hamil- 
tonian of the SC side (a; > 0), which may have a differ- 
ent electronic filling Cn, has to be supplemented by a fi- 
nite (proximity-induced) order parameter A(r) — Aqc*"^, 
where cj) is the superconducting phase. The scattering 
amplitudes of such an interface can be calculated along 
the lines of Ref. [l^l (details are given in Appendix 2). 
Here, we restrict ourselves to subgap transport [based on 
Andreev reflection (AR)] where the quasiparticle excita- 
tion energy e is smaller than the superconducting gap 
Ao. 

In Fig. m we show the behavior of the AR probability 
and the weight of evanescent modes on the left-hand side 
for an incident spin-f electron as a function of M and the 
angle of incidence 9 at the Fermi energy. Similarly to the 
n — n' junction, (see Fig. [2]), the AR exhibits an asym- 
metric behavior with respect to 9, with a maximum for 
large negative M . Again, an evident correlation between 
the presence of evanescent modes and the asymmetry is 
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FIG. 4: (Color online) Angle dependence of the AR probabil- 
ity (left panel) and the amplitude squared of the evanescent 
(hole-like) mode on the HgTe/CdTe QW side (right panel) 
for an incident spin-t electron at a junction between a low-n- 
doped HgTe/CdTe QW and a SC at the Fermi level (e = 0) 
as functions of the mass parameter M, where —C'l = 5 meV, 
— Cr = 1 eV and Ao — 1 meV. The values of the parameters 
A, B, and D are the same as in Fig. [2] Similarly to Fig. [2l a 
clear asymmetry is present for negative values of M. 



observed; compare the left and the right panels of Fig. 0] 
with each other. The case of an incident spin-| electron 
just follows by replacing —9 in Fig. S) 

We remark here that raising the electronic filling of 
the HgTe/CdTe QW at the left-hand side of the inter- 
face alters substantially the asymmetric behavior of the 
scattering amplitudes in both types of junction (normal 
metal and SC). As expected from a junction between 
two metals, asymmetry in such a scattering problem is 
absent. Strong asymmetry is also absent when the linear 
terms are dominant, as in graphene [13, EE] or low-doped 
HgTe/CdTe wells That is the reason that the effect 
has been missed before. 

Hence, the simultaneous existence of propagating and 
evanescent modes leads to the appearance of a type of 
spin Hall effect manifested as a local spin current density 
flowing along the interface. A similar SHE was recently 
predicted at potential steps in the Kane-Mele model in 
graphene [1^ and on the surface of a three-dimensional 
topological insulator ^2^. 

However, the junctions based on HgTe/CdTe QWs dis- 
cussed in our work are, to the best of our knowledge, the 
flrst composite system where a SHE can be observed in 
the absence of structure and bulk inversion asymmetry 
within each subsystem. 



IV. EXPERIMENTAL REALIZATION 

The observation of such an effect along the interface 
is complicated although the use of spin-resolved scan- 
ning tunneling microscopy techniques or optical meth- 
ods might work [25]. Here, we propose an alternative 
way, depicted in Fig. [SJ in order to take advantage of 
the spin- and angle- dependence of the AR probability. 
When a voltage V is applied across the QW/SC contact, 
electrons are injected through the junction from normal 
metal (point-like) reservoirs labeled Rl, R2 in the fig- 
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FIG. 5: (Color online) Schematic of the AR spin splitter: 
Electrons are injected from an unpolarized normal metal 
(point-like) reservoir Rl (and similarly for R2) due to a bias 
voltage V. As discussed before, we can adjust parameters 
such that electrons with spin-J, (blue full line) arrive at the 
junction with an angle of incidence favorable for AR. Hence, 
they are transmitted as Cooper pairs into the SC. This pro- 
cess produces reflected holes with spin- J, collected into the first 
reservoir Rl (blue dashed line). At the same time, this an- 
gle of incidence is less favorable for particle-hole conversion of 
electrons with spin-t (red full line). Thus, these electrons are 
mostly reflected at the interface to a second reservoir (R2), 
resulting in a spin imbalance between Rl and R2. Hence, 
there will be a finite charge current jp from the reservoirs Rl 
and R2 to the SC and there will be a finite spin current 
from Rl to R2. This spin current yields a spin accumulation 
in the two reservoirs in the direction of the bold arrows. 

ure. Then, as discussed in detail above, the presence of 
evanescent modes showing up at the interface acts on 
the incoming wave as a spin splitter whose efficiency de- 
pends on various parameters. A detailed description of 
the working principle of the envisioned device is given in 
the caption of Fig. [SJ 

To quantify the efficiency of the spin splitter, we cal- 
culate the spin conductance, which is a direct measure 
of the imbalance between spin-f and spin- J, carriers in 
reservoirs Rl and R2, induced by the voltage V. It can 
be written as 

Gs(^,^)^^5]a|ree,.r, (4) 

(T = ± 

where |ree,o-P is the probability for an electron with spin 
a and injected from Rl at an angle 9 = arctan(VF/2L) to 
be reflected at the interface. W and L are respectively 
the width of the system and the distance between the 
reservoirs and the interface. The factor 2 accounts for the 
two reflection processes, from Rl to R2 and vice versa. 

We mention here that the contribution to Eq. (jH from 
topologically protected edge states (if the QW has a finite 
width W) is negligible. Hence, our result also applies to 
finite but wide systems where the overlap between edge 
states at opposite edges is weak. In Fig. [51 we addition- 
ally analyze the energy and band gap dependence of the 
performance of the AR spin splitter. It is clearly visi- 
ble that the device works best at small bias voltages and 
that the ideal value of the band gap M changes slightly 
to more negative values as the bias is increased. 
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FIG. 6: (Color online) Bias dependence of the differential 
spin conductance for a QW/SC contact with —Cl = 50 meV 
and —Cr — 1 eV. Gs is plotted for different values of the 
injection energy: eV — 0. (black dots); 0.3 Ao (red crosses); 
0.6 Ao (blue circles), and 0.9 Ao (green triangles). The other 
parameters are Ao = 1 meV, A = 4 eV A, B — —70 eV 
A^, and D — —50 eV A^. The size of the system is taken as 
W = 500 nm and L = 500 nm. It is clearly visible that the 
spin splitter works best at a small bias voltage. 



1. HgTe/CdTe QW/metal interface 

We first consider the interface between a HgTe/CdTe 
QW and a metal modeled as a highly doped HgTe/CdTe 
QW. 

On the left-hand side of the interface, under the choice 
of parameters M < ^^-^^J^^-c^) and < (£; - Cl)^, 
the wave function is a sum of incoming and reflected con- 
tributions, 

\li{x) = *j(a;) + r^'r(x) + f*f;(a;), (5) 

where r (f ) refers to the reflection amplitude of a propa- 
gating (evanescent) wave. 

In the (|£'l t) , \H1 t}) basis and at given energy E and 
transverse wave vector ky, the two-component spinors 
can be written as 



V. CONCLUSION 

To summarize, we have predicted spin-dependent 
transport properties at interfaces between HgTe/CdTe 
QWs and normal metals as well as superconductors, re- 
sulting in an interface spin Hall effect. It has been identi- 
fied that this effect is clearly connected to the coexistence 
of propagating and evanescent modes at the junction. We 
have proposed a setup to functionalize it as an Andreev 
reflection spin splitter. 
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APPENDIX 

In this Appendix, we write explicitly the wave func- 
tions involved in a junction between a HgTe/CdTe QW 
and a normal metal or an s-wave superconductor. Then, 
we define the scattering probabilities plotted in Figs. 2-4. 

We solve the scattering problem for an incoming spin- 
up electron in the conduction band. All the states are 
plane waves and are written as '${x)e'^'^^'^ , where '^{x) is 
a two-component spinor wave function. 



where k = {k^, ky), — \J k\ ~ k^, and = ■\Jky + 
with Kx > 0. In this Appendix, we use the definition 

— — fc|. The momenta ki 2 are defined in Eq. (2). 

The wave function at the metal side reads 

^{x) = mtix) +i-i'i{x), (9) 

with t {t), the transmission amplitude of propagating 
(evanescent) modes and both transmitted wave functions 
are given by 

The momenta and follow from Eq. (2) after the 
substitution C'l Cr. The spin down states are ob- 
tained from the previous ones Eqs. (2)-(4) and Eqs. (6) 
and (7) by substituting k^ — > —k^- 

The scattering amplitudes are obtained using scatter- 
ing matrix theory by matching the left and right wave 
functions and their derivatives at the interface: 

*(a;,y)U^o- = *(a;,?;)U-^o+, 

(9^*(a;,y)U^o- = a^4'(a;, j/)U_^o+- (12) 
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In Sec. III. A, we investigate the behavior of the trans- 
mission probabihty and the amphtude squared of the 
evanescent modes for n-n' and n-p junctions. 

The current conservation imposes the normahzation of 
the transmission probabihty \t\^ by the ratio of incoming 
and transmitted particle current, 



= 1, 



(13) 



where 



OH 

= -2efcicos6'{(i:) + S)(d(k) + M(k))2 

+ A^kl{D - B) - A^{d{\L) + M {]<:))]. (14) 

The transmitted average current has the same form but 
depends on Cr instead of Cl ■ It is possible to solve the 
matching conditions and to find analytical expressions 
for all scattering amplitudes. However, the expressions 
are still too long to be written down here. 

For fixed chemical potentials. Figs. 2 and 3 present the 

dependence of the transmission probability T(6^ M) = 

-t 

\t(6, M) p ~ and of the amplitude squared of evanescent 

modes \i{0, M)p as functions of the incident angle 9 and 
the mass parameter M . 



2. HgTe/CdTe QW/superconductor interface 

We consider now the interface between a HgTe/CdTe 
QW and an s-wave superconductor. We assume the 
superconductivity is induced on the right-hand side of 
the interface by the proximity effect. Thus the effective 
Hamiltonian, Eq. (1), is expanded by particle- hole sym- 
metry and contains the pairing potential A(x, y) = Aoe"'^ 
((j) is the superconducting phase) as off-diagonal matrix 
elements, namely. 



HbcLG 



h{—ihdr 
A*(r) 



A(r) 

'h{—ihdr 



(15) 



with r = {x,y). 

The wave function on the left-hand side of the interface 
is the superposition of electron-like and hole-like quasi- 
particles, 

+ rhe^ruAx)+fhe^fuAx), (16) 

where r^e represents the amplitude of electrons to be 
reflected as electron-like quasiparticles at the interface 
while The refers to the amplitude of an electron-hole con- 
version, namely an Andreev reflection process. The scat- 
tering amplitudes fee and f^e correspond to the reflection 
as electron- or hole-like evanescent modes. 



By solving the Bogoliubov-de Gennes equation 
HBdc"^ = E'^, we obtain the expressions for the four- 
component vectors 



d(k) + M(k) 




(17) 



*r„(a:) = 




(18) 



The momenta k^. 



I d{K) + M{K) 



V 



(19) 



fc,^ follow 



'kf - kl and ^ y- < 

from Eq. (2). The propagating and evanescent hole- 
like wave functions follow from Eqs. (fT7|) and (fT9| by 
the substitution E —E. Note that the electron-like 
excitations exist on the upper two components of the 
wave function and the hole-like excitations on the lower 
two components. 

We focus on transport only below the superconducting 
gap {E < Aq). Then the eigenfunctions on the supercon- 
ducting side contain only evanescent waves 



*(x) = ts+'^s+ix) +ts-1's^ix) 

+ ts>^'ifs>^{x)+ts'_'fs'_{x), 



(20) 



where the spinors read 



*s+(a;) 



/ diks) + M(ks) 

A{kxS - iky) 
[d(k5) + M(ks)]e^('3-^ 
V A{k,S - iky)e'^fi-^^ 



(21) 



"UsAx) 



^s'^{x) 



( d{ks) + M{ks) 
A{kxS - iky) 
Kks)-HM(ks)]e-*('5+^) 
V A{h,s - tky)e-^(^+^^ J 



( d{Ks)+M{Ks) \ 

-iA{-K.xS + ky) 
[d{Ks) + M{Ks)y^^'^^ 
V -*A(-K,s + fcy)e^(^-*) j 



(22) 



(23) 
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/ d{Ks)+M{Ks) 

-iA{-K^s + ky) 
[d{Ks) + M{Ks)]e-'^^+'t'^ 



and we have introduced the fohowing parameters 



kxs 



= arccos(i?/Ao 



1,2 



KxS = 



\/i^S + ky, 

hvp 



(24) 

(25) 
(26) 
(27) 

(28) 



The momenta squared /c|j^ and 



= —k'g2 are defined 
by Eq. (2) after replacing Cl by Cr and setting E to 0. 

The scattering ampUtudes follow from the matching 
conditions of wave functions and their derivatives at the 
interface as written in Eq. p^ . Their analytical expres- 
sions are long and therefore not presented here. 



Also, due to current conservation, the Andreev reflec- 
tion probability is normalized by the particle current of 
incoming electrons and reflected holes, namely. 



\rh, 



■hr 



Jx 



(29) 



where. 



Jx = 



ME) 

= -2ekiCose{{D + B){d(k) 
+ A\l{D - B) - A\d{k) + 



M(k))}. 



(30) 



Due to particle-hole symmetry, the average current car- 
ried by the reflected hole is defined as j^*" = jx{—E). 

Figure 4 shows the asymmetric behavior of 
the Andreev reflection probability Ra{0,M) = 



\rhei0,M)\ 



2 1 jx(-g) I 
Jx(B) 



as a function of the angle of 



incidence and the mass term and its correlation with 
the non-zero amplitude squared of hole-like evanescent 
modes \rhe{d, M)|. 
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